Scaling approach to related disordered stochastic and free-fermion models 
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Motivated by mapping from a stochastic system with spatially random rates, we consider disor- 
dered non-conserving free-fermion systems using a scaling procedure for the equations of motion. 
This approach demonstrates disorder-induced localization acting in competition with the asymmet- 
ric driving. We discuss the resulting implications for the original stochastic system. 
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I. INTRODUCTION 

The role of disorder in stochastic non-equilibrium sys- 
tems is a topic of much recent interest (see e.g., [l|, ^ 0, 
0, [1] and references therein). This is exemplified by stud- 
ies of the asymmetric simple exclusion process (ASEP), 
a one-dimensional lattice gas model exhibiting a phase 
transition for open boundary conditions. The effect of 
quenched particle-disordered hop rates in this system was 
investigated in, e.g., [1, 0, [1| and the case of quenched 
spatial disorder has been clarified by a number of works 
(see, e.g., including via a linearization trans- 

formation and scaling approach jl2i] . 

In this paper we study a related non-equilibrium model 
with applications in catalysis — the partially asymmetric 
exclusion process combined with deposition and evapo- 
ration of dimers (or equivalently pair annihilation and 
creation). Studies of the pure model e.g., [l^, [ij, [l5| 
suggest that, unlike the ASEP, for adsorption and des- 
orption rates which do not vanish in the thermodynamic 
limit, there is no steady-state phase transition even with 
open boundary conditions However, the effect 

of disorder on both steady-state properties (such as the 
density) and dynamics is still of interest. 

Here we demonstrate that, for spatially disordered 
rates obeying a particular condition, one can approach 
this problem via a powerful mapping to free-fermion sys- 
tems. In the pure case, the resulting free-fermion system 
can be treated by Fourier and Bogoliubov transforma- 
tions [H, [l^. The simplest disordered scenario, the case 
of a single defect, is known to produce the stochastic 
analogue of localized modes [19]. One might also expect 
to make progress for the case of dilution, by breaking 
the chain up into finite uniform sections (compare the 
treatment of the diffusion-only problem in [20| ; a similar 
break-up into effective pure segments also occurs in Q). 
In this paper we are chiefly interested in more general 
disordered cases where there is no translational invari- 



ance and therefore one can not utilize the usual Fourier 
and Bogoliubov transformations. However, work by Merz 
and Chalker involving doubling the number of fermion 
operators [2l| , offers the possibility of using a type of a 



localized Bogoliubov transformation followed by numeri- 
cal calculation. Another potential way to make progress 
is to derive and work with the linear dynamic equations 
and this is the approach pursued in the present work. 

Our scaling approach is based on a method developed 
by Pimentel and Stinchcombe [l^] to treat the equation 
of motion of a ID Mattis-transformed Ed wards- Anderson 
Heisenberg spin glass and later applied by the present 
authors to the Cole-Hopf transformed ASEP [12]. Es- 
sentially, one looks at the evolution of parameters after 
a 5 = 2 decimation of the equation of motion. For ex- 
tended states the parameters are found to evolve chaot- 
ically under the scaling whereas for localized states they 
decay exponentially allowing the identification of a local- 
ization length. Application of this method to the present 
problem contributes to understanding of disordered free- 
fermions as well as illustrating a potential way to make 
progress with the disordered stochastic model. 

The plan of the rest of the paper is as follows. In Scc.HIl 
we define our model, introduce the quantum Hamiltonian 
formalism and demonstrate the free-fermion reduction for 
a particular choice of disorder. We also discuss the known 
pure results. In Sec. lIIIi we derive the coupled equations 
of motions for the disordered free-fermion case. Then, in 
Sec. II VI we exploit the fact that these equations are linear 
to apply scaling techniques based on those in [l^, [2^ . 
As in those works, we find that disorder induces local- 
ization effects and in Sec. |V] we tentatively discuss the 
consequences for the original stochastic system. Finally, 
Sec. I VII contains conclusions and suggestions for future 
work. 



II. MODEL AND MAPPINGS 
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Let us start by considering the pure partially asymmet- 
ric exclusion process with dimer evaporation and depo- 
sition, shown schematically in Fig. [TJ Subject to the ex- 
clusion constraints, particles hop to the right (left) with 
rate p (q), and pairs of particles are deposited (evapo- 
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FIG. 1: Partially asymmetric exclusion process with dimer 
evaporation and deposition (pure model). Filled circles indi- 
cate particles; open circles denote vacancies. 



rated) with rate e (e'). This provides a simple model 
for adsorption-desorption processes and catalytic surface 
reactions. 

In the special case p + q ~ e + e', the system can 
be mapped via a Jordan- Wigner transformation [23j and 
a Bogoliubov-Valatin transformation [23, ^2^ to a free- 
fermion problem and hence solved exactly ^131 3 
The results are of more general significance since, by du- 
ality transformations and spin rotations, the quantum 
Hamiltonian can also be mapped to other stochastic sys- 
tems 15, 271, including some which are experimentally 
realizable ^] . A subcase is p = q = e = e' which maps 
to the Ising model and so is well understood p^ . 

The aim of the present work is to study the spatially- 
disordered version of this model,with rates pi, qi, Si and 
e'l for processes on the bond between sites I and I + 1. 
In this section we demonstrate that the mapping to free- 
fermions is still valid in the case where pi + qi = £i + £[ 
for each I. 

The generalized pair process parameterized by pi, qi, ei 
and e'l can be conveniently represented using the quan- 
tum Hamiltonian formalism [2^. In this approach one 
defines a probability vector |P) = '^'^^^ I") 

the basis vector associated with the particle configura- 
tion n = (ni, 712, . . . , n^) and P„ the probability mea- 
sure on the set of all such configurations. \P) obeys the 
normalization condition {s\P) — 1 where (s| = X]n("'l 
and (ri|n') — bn,n'- Within this formalism the master 
equation for the time evolution resembles a Schrodinger 
equation 



-H\P(t)) 



with 



- \{vi +qi+£i+ e'l) + ji-pi +qi+£i- e'i)al 



+ ^{pi-m+£i- e'Mt+i +'^{pi+qi~£i- e'lhfcrf^i] ■ 



(s\H = 0. 



(3) 



Note that the general four-parameter model defined 
by ([2]) contains many other important models as sub- 
cases, e.g., the ASEP {qi — ei = e'l — 0), symmetric 
diffusion {pi — qi, ei — e[ — 0), etc. Further progress 
depends on what boundary terms are imposed and how 
tractable the resultant quantum spin model is. In partic- 
ular, a Jordan- Wigner transformation to fermion opera- 
tors cj, ci gives both an easily treatable quadratic part 
(arising from and terms) and a more difficult 

quartic part (from the a^a^ terms). The quartic terms 
will clearly be zero for the special case 



pi + qi = ei + £; 



(4) 



and this gives the disordered free-fermion model which 
will be explored in more detail in the remainder of this 
paper. 

For this particular case, we now outline the details of 
the transformation to fermion operators (compare Gryn- 
berg et al. [18] for the pure model). Writing cr^ in terms 
of cr+fT^ and imposing periodic boundary conditions, 
Eq. ^ becomes 

I 

+ [Pi-i -Pi+£i- e'i_^)a+ai - (5) 

Then applying the usual Jordan- Wigner transforma- 
tion [H 



cr+ = c|exp I i7i-^c]( 

3<l 



CFi = CXp I -ZTT^cjcj I Q, 
3<l 



(6) 



(7) 



(1) yields the Hamiltonian in terms of fermion operators 
Hff = ^[p/Qc|+i - qicjci+i - eiicjcl^^ - 1) 



I 

+ e'icici+i - (p/_i -pi+ei - e'i_i)clci]. (8) 

The study of the disordered free-fermion Hamiltonian ([8|) 
is the central aim of this paper; we now digress to sum- 
marize the known pure results. 

In the pure case one can exploit translational in- 
variance via a Fourier transformation and then use a 



(2) 



Here the "spin-flip" operators cr^ correspond to particle 
creation and annihilation processes while the projection 
operators az are needed to account for the probability 
that configurations do not change. Conservation of prob- 



^ If the total number of fermions, ^[O'^'^^ = cjci, is even then 
the fermion operators must obey cycUc boundary conditions; in 
the odd subspace, anti-cycUc boundary conditions are required. 
See [H. 
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Bogoliubov-type similarity transformation [2J, [25| (well 
defined only for e and e' non-zero) to diagonalize the 
Hamiltonian. The details are given in [isl. Il8l|: here we 
simply quote the resulting Hamiltonian: 



(9) 



with spectrum 

Xk = e + e' + {e — e') cos k + i{p — q) sin k. (10) 

The imaginary part of Xk implies ballistic motion while 
the real part indicates that excitations decay with time 
constant [e + e' + (e — e')cosfc]~^. For e and e' non- 
zero, the spectrum is gapped giving exponentially fast 
kinetics. By mapping back from this free-fermion model, 
one can obtain the steady-state density profile for the 
original pure stochastic problem. The result 



Pi 



(11) 



is in agreement with mean-field calculation (one can show 
that the system is spatially uncorrelated). Other works 
have calculated dynamic correlation functions [l^, [3], 
shock evolution [30,], persistence probabihty, etc. Note 
that non-steady-state properties of the stochastic system 
can be understood by mapping to the Glauber model [sTIl 
and treating the asy mmetric hopping as biased diffusion 
of domain walls p32l |. 

In the disordered case where there is no translational 
invariance a different approach is required — in the next 
section, we demonstrate how to derive the equations of 
motion starting from the free-fermion Hamiltonian ([8|l. 



III. EQUATIONS OF MOTION 

Let us first consider the expectation value of a general 
observable X, i.e.. 



(X) (i) = (,s|X|P) 



(12) 



where we use angular brackets to denote this stochas- 
tic average over histories. Using the "Schrddinger equa- 
tion" ([T]) one sees that the time evolution of {X) is given 
by 



d{X) 
dt 



= {[H.X]) 



(13) 



with [H^X] the usual commutator. 

For the free-fermion Hamiltonian ([5]) we can use the 
fermion anti-commutation relations, together with the 
identity [A, BC] = {A, B}C - B{A, C}, to calculate the 
commutators [cj , H] and [cj , H] and hence obtain coupled 
dynamic equations for the operators (in the Heisenberg 



representation) : 
dcj 



+ iP]-! -Pj+- 



(14) 



dt 



{Pj- 



Pj 



(15) 



Note that due to the non-Hermicity of B these equations 
are not simply Hermitian conjugates of each other. The 
equations are linear thus raisin g th e possibility of using 
scaling techniques like those in [12, 124I — this will be the 
subject of the next section. 

As a preliminary to this scaling approach, we first as- 
sume that Cj and cj can be written as superpositions of 
operators with exponential time dependence: 



c,=^r,(c.)e"*, 

UJ 



(16) 
(17) 



Substituting into Eqs. and (fT5|) and equating com- 
ponents with the same time dependence yields 

(Ej - uj)Vj = -pj_irj_i - ^jFj+i -I- £j-i7j-i - ej7j+i, 

(18) 

{Ej -f w)7j = -Pj7j+i - gj-i7j-i + ej-irj-i - Sj-rj+i, 

(19) 

with = (pj_i - -h Ej - e^_i). 

Note that for the pure case we can carry out a Fourier 
transformation by writing Y j(Lo) = e^^^^Tk, Jj{<jj) = 
C'^^-'^k- As expected, we then find that for each Fourier 
component, Eqs. (fT8|) and (fT9|) have solutions with cj = 
Afe or w = — A_fe, where is the pure spectrum ITUl) . In 



other words the operators Cj and can be written as 



'ikj 



(Ffee 



where 



e' k 
Ik = iFfc— cot -, 
£ 2 



ik = -*rfctan-. 



(20) 
(21) 

(22) 
(23) 



The Bogoliubov angle can be obtained from the rela- 
tions dm and P5|) : in the pure case this approach thus 
provides an alternative route to direct transformation of 
the Hamiltonian ([5]). 

As an aside, we remark that Eqs. (fT5|) and can 
also be cast in the form of a transfer matrix mapping. 
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Analysing products of disordered transfer matrices then 
provides a way to investigate possible localization effects. 
This is a particularly attractive approach for simple mod- 
els such as binary disorder. In the next section we de- 
velop instead, a general numerical scaling approach for 
arbitrary distributions of disorder. 



IV. SCALING 

A. Procedure 

In the general disordered case it is not obvious how to 
solve Eqs. (fT5|) and (fTO)) explicitly but their form is rem- 
iniscent of an equation which appears in the disordered 
ASEP (in a mean-field description) after a linearizing 
(Cole-Hopf) transformation so we can develop a scaling 
method similar to the one used in that case [I4I (compa.re 
also the original work of Pimentel and Stinchcombe [22|). 
For full generality we must allow each coefficient to scale 
differently so we perform a, b — 2 dilation of the system 
using decimation on the following equations 



-^ij+i7j+i, (24) 



(25) 



where we have extended the notation for clarity. For 
example, Pjj-i denotes the coefficient of Fj-i in the ex- 
pression for Tj. Even in the pure case, this may scale 
differently from Pjj+i which is the coefficcnt of 7jj+i 
in the expression for 7^. In the rescaled equation for Tj 
(7j) the coefficient of rj_2 (7^+2) isPjj-2 (Pjj+2)- The 
resulting 10-parameter scaling equations are rather in- 
volved and are displayed for convenience in Appendix [Xl 
We now test this scaling procedure on the pure model 
(Sec. IIVB|) before applying it to the disordered problem 
fSec UVCj) . 



B. Pure model 

We here demonstrate how the scaling of the equations 
of motion and P5|) reflects the known spectrum Afc 
of the pure model ^TU\i . For definiteness we consider the 
case e > e', p > q throughout; the extension to other 
cases should be straightforward, since they are related 
by duality or space reflection. 



1. Non-biased case, p ~ q 

For the pure case, if one starts with p = p = q = q, 
e = e, e' = e' , and by implication E = E, then the rela- 
tionships p — q, p — q, e = e and e' = e' are preserved 
under scaling so one only needs to scale six independent 




FIG. 2: Evolution of parameters in Eqs. p4|l and (|25|) under 
6 = 2 decimation. Pure case with uj = 0.41 and starting 
parameter values p — p — q — q = 0.5, e = e — 0.8, e' = 
e' = 0.2, E — E = 0.6. Graph shows absolute magnitude of 
parameters against number of iterations /; lines are provided 
as an aid to the eye. Chaotic behaviour is typical of scaling 
in the allowed band. 



parameters. We find two distinct types of scaling be- 
haviour. 

For 2e' < ijj < 2e the scaling parameters all evolve 
chaotically under scaling (see Fig. [2]). In contrast for 
w < 2e' or w > 2£, the "energies" E and E tend to con- 
stant values whereas the rescaled rates r (i.e, p, g, etc.) 
tend rapidly to zero inside an exponentially decaying en- 
velope (see Fig. [3]), i.e., r{l) ~ /(O^ '^^ where ^ = 2^ is 
the distance between sites (/ is the number of iterations) 
and ^ is a localization length. We also see this localized 
behaviour for all non-real uj. 

Just as in the scaling of the Cole-Hopf-transformed 
mean- field ASEP [12] the localization length can be mea- 
sured numerically from the exponential decay of the pa- 
rameters. In Fig. m we plot the localization length ob- 
tained from the scaling of p, as a function of real fre- 
quency UJ (localization lengths obtained from the other 
rates are found to be identical). These numerical results 
can be simply explained by reference to the pure spec- 
trum for the symmetric case. The "allowed" band (real 
k) corresponds to the region 2e' < uj < 2e\ here states 
are extended and the localization length is infinite. Fre- 
quencies outside this band correspond to a complex wave- 
vector k+iK, i.e., localized states with localization length 
^ — Indeed, from uj — e + e' + {s — s') cos(fc -I- in) we 
can obtain an analytical expression for the localization 
length on the real-o; axis (fc = or A: = tt): 



^ = I arccosh 



- e — UJ 



(26) 



This expression is shown as a dashed line in Fig. [?] and 
agreement with the numerics is excellent. 
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FIG. 3: Evolution of parameters in Eqs. (|24p and (|25|l under 
b = 2 decimation. Pure case with uo = 0.39 and starting 
parameter values p — p — q — q = 0.5, e — e = 0.8, e' — 
e' = 0.2, E = E = 0.6. Graph shows absolute magnitude of 
parameters against number of iterations I; lines are provided 
as an aid to the eye. Rate parameters decay exponentially as 
a function of 2^, as expected outside the allowed band. 
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FIG. 4: Localization length ^ versus (real) frequency to for 
pure case with p = q = 0.5, e = 0.8, e' = 0.2. ^ is effectively 
infinite in allowed band 0.4 < u < 1.6. Crosses indicate data 
from numerical scaling of rate p; dashed line is analytical 
prediction of Eq. pSjl . 



2. Biased case, p ^ q 

In the pure biased case the "allowed" values of oj are 
given by Xk of Eq. pO|) with k real — an ellipse in complex 
Afe-space, see Fig. [5] For these values of lu we again find, 
as expected, that the scaling parameters all evolve chaot- 
ically under scaling corresponding to extended states. 
Outside this allowed ellipse we see localized states once 
more. 

Inside the ellipse we see a new type of behaviour — E 



ImAfe 



localized 




ReXk 



FIG. 5: ImAfe versus ReA^ for pure model (case p > q, e > 
e'). Solid ellipse represents allowed spectrum given by (|10|l . 
Inside the ellipse equations of motion scale to ballistic limit 
p/q — > oo; outside of ellipse is localized region. 




FIG. 6: Evolution of parameters in Eqs. (I24p and (I25p under 
6 = 2 decimation (first few iterations). Pure case with uj = 
0.41 and starting parameter values p — p = 0.55, q — q — 
0.45, £ = e = 0.8, e' = e' = 0.2, E = E = 0.6. Graph 
shows absolute magnitude of parameters against number of 
iterations I; lines are provided as an aid to the eye. Rates p, 
q, e, e' scale exponentially to infinity while p, q, e, e' decay 
exponentially and "energies" E, E tend to constant values. 
This is typical of ballistic behaviour inside the allowed ellipse 
(for case p > q). 



and E again tend to constant values, some of the rates 
scale exponentially to zero, whereas others grow expo- 
nentially to infinity (see Fig. [6]) . We argue that in this 
region the system is dominated by the ballistic motion 
and scales to the totally asymmetric case. So if we start 
with p > q, then p scales to infinity while q decays to 
zero. 

In Figs. [7] and [8] we plot the localization lengths for 
p and q along the principal axes of the w-space ellipse 
(again for the case p > q). As is evident from the dis- 
cussion above, the localization length for q is small and 
positive in all regions — q is an irrelevant parameter under 
scaling. The localization length for p is positive in the lo- 
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FIG. 7: Numerical localization lengths ^p, as a function 
of real frequency ujrc for pure case with p = 0.55, q — 0.45, 
£ = 0.8, e' — 0.2. Note difference in scaling behaviour for p 
and q. Lines are provided as an aid to the eye. 
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FIG. 9: Evolution of parameters in Eqs. pljl and (f25|l under 
b — 2 decimation. Disordered case with u — 1.41 and starting 
parameter values pj.j-i = Pj-i.j = qj-i.j = =0.5, 
Sj.i-i = ^i-iJ = ^1 = = '^-z, E = E = Pj,j-i - 

Pj+i,j — £j+i,j — s'j.j-i where z is a random variable drawn 
from a uniform distribution between 0.6 and 1.0. Graph shows 
absolute magnitude of parameters at site jo (averaged over 
100 realizations of disorder) against number of iterations I; 
lines are provided as an aid to the eye. Rate parameters decay 
exponentially as a function of 2^ showing a clear localization 
effect. 



the case where the disorder is in e and e' with the rates 
chosen to maintain the free-fermion condition for each 
bond (specificaUy e'j = 1 — Sj, pj = p and qj = q, where 
p+q — I)- Once again, we consider in turn the symmetric 
(p = q) f-iid asymmetric (p ^ q) cases. 



-150 ' ' ' ^ ' ^ ' ' ' 

-0.4 -0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 0.4 

FIG. 8: Numerical localization lengths ^p, as a function 
of complex frequency uj = 1 -\- iuJim for pure case with p = 
0.55, q — 0.45, e — 0.8, e' — 0.2. Note difference in scaling 
behaviour for p and q. Lines are provided as an aid to the 
eye. 



calized region outside the eUipse, infinite on the allowed 
ellipse, and negative in the ballistic region inside (refer 
again to Fig. O. The value of this relevant localization 
length can again be understood analytically by consider- 
ing a complex wavevector in the dispersion relation. 



C. Disordered model 

Armed with this understanding of the pure case, we 
move on to look at the effect of disorder on the free- 
fermion system. There are obviously various different 
possible ways to add disorder; here we mainly investigate 



1. Non-biased case, p = q 

In Fig. [5] we show an example of the scaling behaviour 
for the case with pj = qj = 0.5, Sj drawn from a uniform 
distribution and e'j — 1 — Sj . The value of lo corresponds 
to being in the allowed band for the pure case (compare 
Fig. (2) . We find that the rates evolve chaotically towards 
zero within an exponential envelope and so we can define 
a localization length via r{l) ~ f{l)e~^^^. In Fig. [TUl we 
plot the localization length for the decay of p, averaged 
over realizations of disorder. Localization is clearly seen 
for all values of lu in the pure band. 



2. Biased case, p ^ q 

In Fig. [11] we plot the localization lengths for p and 
q along the real frequency axis for a sample disordered 
case, viz., pj — 0.55, qj — 0.45, Sj drawn from a uniform 
distribution and e'j — I — Sj. The pure results are also 
shown for reference (compare also Fig. [7]). In this case we 
see that disorder reduces the size of the allowed ellipse 
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FIG. 10: Localization length for p scaling versus (real) fre- 
quency for disordered case with p = q = 0.5, Sj = z, e'j = 1 — z 
where z is drawn from a uniform distribution between 0.6 and 
1.0. Pure results outside the band are also shown for compar- 
ison (dashed). In the disordered case, the localization length 
is finite everywhere except uj = 1. 



FIG. 12: Numerical localization length as a function of real 
frequency ljro for disordered cases with Sj = z, e'j = 1 — z 
where z is drawn from a uniform distribution between 0.6 
and 1.0. Results are shown for a range of values of p and q as 
indicated by the legend. Lines are provided as an aid to the 
eye. 
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-100 ' ' * ^ • ' • ^ * • ' 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

FIG. 11: Numerical localization lengths ^p, as a function of 
real frequency ujrc for disordered case with p = 0.55, q — 0.45, 
Ej = z, e'j = 1 — z where z is drawn from a uniform distribution 
between 0.6 and 1.0. Lines are provided as an aid to the eye. 
Pure results are shown for comparison. Note the reduction of 
the ballistic region in the disordered case. 



along the real-w axis, i.e., it acts to increase the localized 
region at the expense of the ballistic region. The results 
for CO — 1 + iuji^ are essentially the same as in the pure 
case (Fig. H). 

This competition between disorder-induced localiza- 
tion and driving is analogous to that arising in another 
non-Hermitian Hamiltonian studied by Hatano and Nel- 
son (33I I . Their work was motivated by the mapping be- 
tween flux lines in a (ci-l-l)-dimensional superconductor 
and d-dimensional bosons. A random potential in the 



bosonic problem arises from the study of columnar de- 
fects in the superconductor. 

For all cases with disorder only in s and e' , we find that 
the localization length is altered along the real-o; axis but 
unchanged along the axis in the imaginary direction.^ 
However, the change in the position of the allowed band 
(where the localization length is infinite) depends on the 
strength of the driving p ~ q, see Fig. [T^l We shall now 
consider some tentative analytical ideas to quantify these 
disorder effects. 

Recall that the pure spectrum is given by 

^ ba — gq cos k + i{p ~ q) sin k, (27) 

with bo = e' + e and gq = e' — e. We seek to modify 
this to represent the effective spectrum of the disordered 
model. For the case of disordered e and e' (with £j +£j- = 
p + q ^ 1) then a naive suggestion is 

= bo — a cos k + i{p — q) sin k (28) 

where a can, in principle, have both real and imaginary 
parts, i.e., 

a^d + i2S (29) 

with a, S real. The imaginary part of a causes a rotation 
of the allowed ellipse (see Fig. [121); if 5 7^ oq, then there 
is also stretching or compression. In agreement with the 
numerical localization length data, this form leaves the 



^ Obviously disorder in p and q affects the imaginary-o; direction 
instead. 
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FIG. 13: Schematic of effective disordered dispersion relation 
in w-space for Xk = bo + (ao + i2S) cos k + i{p — q) sin k (solid 
line). Pure case Afe = 60 + ao cos k + i{p — q) sin k is shown 
for comparison (dashed line). Inside the ellipse equations of 
motion scale to ballistic limit p/q — > 00; outside of ellipse is 
localized region. 



centre of the ellipse and the intersection with the lo = 
bo axis unchanged, It also explains why any amount of 
disorder causes localization in the p — q case (where the 
ellipse is a line along the real axis) for all values of real 
frequency except uj — 1. 

There is no a priori reason why the imaginary part of 
a should have a positive sign so we naturally expect the 
possibility of both ±2(5. Furthermore, we anticipate that 
d is related to the mean of e' — e and d is proportional to 
the width of the distribution. In fact, comparisons of the 
obtained localization lengths (e.g., see Fig. [72]) with the 
prediction from the dispersion relation suggest that 
the scaling results can not be described by unique values 
of a, S. It appears that what we measure in the numerical 
scaling procedure, is the average of 1/^ for some range 
of values. Further work is needed to check and quantify 
this hypothesis. 



V. DISORDER EFFECTS IN STOCHASTIC 
SYSTEM 

In this section we attempt to infer disorder effects 
in the original stochastic hopping process by mapping- 
back the results for the equivalent quantum free-fermion 
model. In particular, we are interested in disorder- 
induced changes to the pure steady-state density pT|) . 

We argue from the previous section that one can 
crudely characterize the effects of disorder in the 
evaporation-deposition rates by the replacements e' —> 
e' ± iS, e ^ £ ^ i6 (oq ^ oq ± i25) where 6 is related 
to the width of the disorder distribution. This gives two 
different effective Bogoliubov angles (possibly connected 
to the doubling of fermionic degrees of freedom in the 
work of Merz and Chalker [IH). Stochastic observables 
are presumably related to the average of these two pos- 
sibilities but it is not yet clear what is the appropriate 
function of e, e' to average over. Motivated by the form 



0.6690 
0.6685 
0.6680 
0.6675 
0.6670 



0.6665 



p = 0.5, g = 0.6 
p = 0.55, q = 0.45 
p = 0.7, q = 0.3 
p = 1.0, q = 0.0 
PO + 



0.00 



0.05 



0.10 
W 



0.15 



0.20 



FIG. 14: Simulation results showing density dependence on 
width of e, e' distribution for range of p, q values (see legend) 
in a system of size 1000. Sj is drawn from a uniform distribu- 
tion between 0.8 — w and 0.8 -I- w; e'j given by 1 — Sj. Dashed 
line is a squared fit given by p = po + bw^ with b = 0.05. 



of the Bogoliubov transformation [cf. Eqs. (f22|) and (j23|) ] 
and by comparisons with simulation, we suggest that the 
ratio e'/e is the important quantity and hence define 



e' -id 



— iS e + i6 
This then yields a disorder-averaged density 

_ 1 

Pi = 



1 



(30) 



(31) 



which is increased from the pure case by an amount pro- 
portional to (5^ (for small 5) . Simulations on the stochas- 
tic model (see Fig. [T^ confirm a density shift propor- 
tional to the square of the width of the distribution, i.e., 
the localization in the disordered free-fermion model does 
appear to map back to a density shift in the equivalent 
disordered stochastic system. Further work is still needed 
to rigorously establish this picture and to relate the value 
of 5 both to the details of the disorder distribution and 
to the numerical localization data. 

Similarly, from knowledge of the effective dispersion 
relation for the disordered case, it should be possible to 
adapt the pure treatment of [11] for both steady-state dy- 
namic correlation functions and non-steady-state proper- 
ties. These quantities typically contain time-dependences 
of the form e^^'^'*'''* so any decrease in the effective spec- 
tral gap would be expected to lead to a slowing-down of 
the dynamics (as was found for the ASEP in p^). 



VI. DISCUSSION 

In this paper we have studied the effect of disor- 
der on a model consisting of the partially asymmet- 
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ric exclusion process combined with dimer evaporation- 
deposition. We presented an exact mapping to an equiv- 
alent quantum system and then employed a numerical 
scaling technique on the quantum equations of motion. 
This provided a clear demonstration of disorder-induced 
localization acting in competition with the asymmet- 
ric driving. We discussed tentatively how these local- 
ization effects are related to the observed steady-state 
density shift in the original stochastic system; exact de- 
tails of the inverse transformation remain to be clari- 
fied. The stochastic system is relatively robust to dis- 
order since particle pairs can be evaporated/deposited 
along the whole length of the chain. 

There is much scope for further work on this disordered 
model. In particular, a more rigorous development of 
the mapping from the disordered quantum model back 
to the stochastic system is needed. This would permit 
the translation of numerical scaling results for different 
forms of disorder (including, for example, the case where 
all rates are disordered and Ej + e'j — pj + qj is position- 
dependent) into quantitative statements about the effect 
on density and correlation functions. As indicated at the 
end of Sec.|V]we should also like to extend the discussion 
to dynamics. 

Another open question relates to the effect of disor- 
der in the open-boundary case when the evaporation and 



deposition rates scale with system size, i.e., they are pro- 
portional to 1/L where L is the number of sites. In this 
case the interplay between boundary and bulk effects is 
expected to give a steady-state phase transition in the 
pure model — compare recent work for a similar dimer 
model |34l and earlier results for the case with adsorp- 
tion/desorption of monomers fsij. The addition of 
disorder to the latter pure monomer model was recently 
considered in the context of minimal models for intercel- 
lular transport [37| . 
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APPENDIX A: 



DISORDERED FREE-FERMION 
SCALING 



Here we present scaling relations for the coupled free- 
fermion equations of motion ([24| and ([25|) : 



J 



(Al) 
(A2) 



A & = 2 decimation of these equations leads to 



{Ej +Uj)-fj = ~p'j.j+2lj+2 - qj.J-2l3-2 + {^'yj.j-2^3-2 - i£%,j+2^3+2, 



(A3) 

(A4) 



with rescaled "energy" parameters 



E, = E, - 



^3 - ^3 - 



Pjd--Llj-K 



E. 



E, 



'3+1 



E.-i + uj 



-l£ 



E 



3+1 J 



'3+1 



UJ 



Q3.3-iPj-i.3 P3,3 + i^j+l,3 , ^3,3-1^3-1,3 , £' 3,3 + ^^ 3 + ^,3 



E. 



E. 



E. 



i-1 



E^ 



'3+1 



(A5) 
(A6) 
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and the rescaled rates 



where 



,j+2 



(£')j-,j+2 = 



^ Aj Bj 



{E'.-u>){E'.+u>) 



1 - 



; + lPj + l,j + 2 \ _ Aj / £'j,j + iej + l,j + 2 _|_ P3,j + lPj + l,j + 2 



1 - 



'^'j,j + lgj + l,j + 2 I P3,3 + lP3 + l,3 + 2 \ _ ^3 ( 93,3 + 1^3 + 1,3 + 2 , Sj, j + + j + 2 
^3 + 1-'^ Ej + i+W 7 ^^■+'^ V ^3 + 1-^ Ei + i+U! 



1 - 



_Ej_l-0 



(B^.-w)(i5^.+a;) 

J-igJ-i.j-a A _|_ Bj / Pj,j-i£j- 



-1,3-2 



I £3.J-lgi-l,3-2 A 



1 - 



'j,j-iPj-i.j-2 I jj.j-iej-i 



(£;j-c^)(£;j+(^) 



;iS-i,j-2 A I Bj f Pj,j-iPj-i,j-2 I ej,j-iej-i,j-2 \ 
]^ A 3 Bj 



{E'.-u>){E'.+u3) 



'■3 3 ' ' 

'3,3 + 113 + 1,3+2 I §3,3 + 1^' 3 + l,3 + 2 \ _ Bj f 9^,3 + 193 + 1,3 + 2 i £3,3 + 1 ^'3 + 1,3 + 2 \ 

-^3+1-" gj+i+g- ) V -Ej+i-'^ gj+l+UJ / 
AjBj 



(E'^-U3){E'.+U3) 



E3-1 - 



ij + l + LO 



^ ' E'j_i-w E'^'+i-w Ej^i+u Ei 



5j+i + tj 



Note the difference in meaning between e' (the evaporation rate) and (eY (the scaled deposition rate). 



(A7) 
(A8) 
(A9) 
(AlO) 
(All) 
(A12) 
(A13) 
(A14) 



(A15) 
(A16) 
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